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Kiihn and Mazzeo Reply Van Enter, Kiilske and 
Maes (EKM) in their comment Q to our letter Pj point 
out - correctly - that joint probability measures on 
the product of spin- and disorder-space in systems with 
quenched randomness are non-Gibbsian. They conclude 
from this - incorrectly in view of the connotations carried 
by such a characterization - that our approach which pro- 
vides approximate descriptions of such systems in terms 
of a net of Gibbs measures suffers from being ill defined. 

To wit, it has been shown J3| that the approximation 
scheme used in provides a variational family of in- 
creasing exact lower bounds for the free energy of systems 
whith quenched randomness in terms of Gibbs measures 
on the joint spin- and disorder-space. These measures 
agree on an increasing set of moments with the mea- 
sure characterizing the quenched disorder. Moreover, at 
each level of the approximation studied in ||,0J the vari- 
ational problem is convex (also in the thermodynamic 
limit), hence its solution unique. In this sense, our ap- 
proximation scheme is perfectly well defined, and it is 
legitimate to investigate its efficiency in analyzing un- 
solved problems in the physics of systems with quenched 
randomness. 

The joint measure related to the 2D spin-diluted 
Ising model studied in is Gibbsian in the high- 

temperature phase, while only weakly Gibbsian in the 
ferromagnetic low-temperature phases (for definitions see 
the literature cited in |Q). Given the representation of 
the approximate measures studied in , it is perfectly 
conceivable that they could approach the joint measure 
of the quenched system, if one were able to carry the 
approximative scheme to its end. While a formal proof 
of this statement may be difficult, it nonetheless indi- 
cates that differences between the limiting measure of 
our scheme and the weakly Gibbsian measure it attempts 
to approximate in the low temperature phases might be 
rather subtle. 

Rather than formal convergence properties, which ulti- 
mately require clarification, our main issue from a prag- 
matic point of view is the efficiency of the approxima- 
tions in describing physical properties. As is quite com- 
mon, we are only able to work through the first few 
steps of the scheme, hence the pragmatic question: Is the 
physics we are thereby missing essential or not? As far as 
one can presently see, it is not. Our findings concerning 
critical phenomena in 0,0] are fully in line with currently 
available empirical evidence coming from other conven- 
tional approaches, as well as with known exact results. 
The fact that critical behavior is, nevertheless, still un- 
der debate for the 2D spin-diluted Ising model is related 
to the difficulty in disentangling the finite size signatures 
of the conflicting scenarios from the numerical data, irre- 
spective of what kind of (efficient) approximation is being 
used to obtain them H . 

The physics we are likely to miss generally is that as- 
sociated with Griffiths' singularities 0]. We have been 


able to prove this for the II? system ||: the simplest 
approximation in this case (system (a) of [||) is exact in 
zero field h = 0, and it also very accurately describes 
critical behavior and scaling in the vicinity of the mul- 
ticritical point T = h = 0, p = I, but fails to exhibit 
Griffiths singularities. The absence of these very weak 
essential singularities, while certainly a drawback, docs 
in the light of available evidence not appear to be a very 
serious one for the problem studied in H0J. 

The important message of EKM is that formal conver- 
gence issues of our approach have indeed not yet been 
dealt with in a satisfactory way. A number of such issues 
occurring in situations where non-Gibbsian measures do 
play a potentially worrying role have been studied before, 
mainly in the renormalization group (RG) context (see 
in particular || , and other references in [Q). 

However, our case is also different from the RG set- 
ting in essential aspects, the perhaps deepest one being 
that the disorder potentials used to express approximat- 
ing measures are invariant under local spin-reversal, a 
feature which may well remove many of the pathology 
generating mechanisms occuring in the RG context. 

In summary, there is clearly nothing automatic about 
the success of our approach, nor, however, about its po- 
tential failure due to some formal similarities with the 
RG situation. Independent evidence must in any new 
case be checked. We have done no less in with 
clearly encouraging results. 

We are indebted to Charles Newman for very useful 
discussions. G.M. acknowledges support from the Foun- 
dation Blanccflor Boncompagni-Ludovisi. 
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